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In this paper we investigate quantum metastability of a particle trapped in between an infinite wall 
and a square barrier, with either a time-periodically oscillating barrier (Model A) or bottom of the 
well (Model B). Based on the Floquet theory, we derive in each case an equation which determines 
the stability of the metastable system. We study the influence on the stability of two Floquet states 
when their Floquet energies (real part) encounter a direct or an avoided crossing at resonance. The 
effect of the amplitude of oscillation on the nature of crossing of Floquet energies is also discussed. 
It is found that by adiabatically changing the frequency and amplitude of the oscillation field, one 
' can manipulate the stability of states in the well. By means of a discrete transform, the two models 

are shown to have exactly the same Floquet energy spectrum at the same oscillating amplitude and 
frequency. The equivalence of the models is also demonstrated by means of the principle of gauge 
invariance. 
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I. INTRODUCTION 

> 
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Recent developments in powerful lasers have induced great theoretical and experimental interests in quantum 
systems under the influence of periodic driving fields. For example, an active topic in the field of atomic physics is 
the dynamical stabilization of atoms by intense high-frequency laser fields 0, 0, 0, 3 • Scattering of particles by a 
time-periodic potential is frequently used as a model for the photon-assisted tunneling and other quantum transport 
problems 0, S S S 0- The interesting phenomena of coherent destruction of tunneling in a periodically driven 
two- level systems has also been discussed extensively in the last decade [IollrT| . 

Although diverse phenomena related to time-periodic potentials have been widely investigated, in our opinion the 
Oh, interesting issue of quantum decay through tunneling of a metastable system trapped in a time-dependent potential 
' should receive more attention than it has thus far. Metastable systems which decay through tunneling exist and are of 
interest in many areas of studies, the most famous of which being the ev-decay of atoms. In many realistic situations 
of decay processes, the metastable systems have to be treated as being influenced by time-dependent potentials. 
q-( For example, in some versions of the inflationary models of the early universe fl^ . inflation is governed by a Higgs 
• • i field trapped in a metastable state. Inflation ends when the metastable state decays to the true ground state of the 
universe. During inflation the universe expands exponentially. It is thus obvious that the metastable state of the 
Higgs field is trapped in a rapidly varying potential. Unfortunately, in these inflationary models the true forms of the 
time-dependent potentials are determined by the dynamics of the system itself. This inherent difficulty makes the 
complete analysis of time-dependent decay of the Higgs field rather difficult. It is therefore desirable to gain some 
insights first by studying metastability in various types of time-dependent potentials in simple quantum-mechanical 
models. 

In our previous work |l3j | we have proposed a simple metastable system with a moving potential which has height 
and width scaled in a specific way introduced by Berry and Klein [14| . In that model we found that a small but finite 
nondecay probability could persist at large time limit for an expanding potential. Another class of time-dependent 
metastable system had been discussed by Fisher [l^, in which he considered a generic class of static metastable well 
subjected to a weak time-periodic force. As the strength of the time-periodic force was taken to be small, it could 
be treated as a small time-dependent perturbation term. Fisher modified the standard (time-independent) WKBJ 
approach to include a weak time-dependent perturbation. He found that the weak time-periodic force always enhances 
the time average of the decay rate of the system. The potential considered by Fisher has a number of oscillator-like 
levels near its minimum. The opposite situation where only two levels are present was considered by Sokolovski 
[l6|. who studied the effect of a small AC field mixing two levels in the well on the tunneling rate in a semiclassical 
framework. 

Recently, we have considered the problem of how the metastability of a quantum system used to model the a-decay 
was influenced by an oscillating barrier [l7|. The Floquet formalism was adopted in our analysis, so that the periodic 
field needs not be treated perturbatively, and the number of states in the potential is not restricted. Our results 
show that an oscillating potential barrier generally makes a metastable system decay faster. However, the existence 
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of avoided crossings of metastable states can switch a less stable state to a more stable one, and vice versa. It is also 
found that increasing the amplitude of the oscillating field may change a direct crossing of states into an avoided one. 
Hence, by a combination of adiabatic changes of the frequency and the amplitude of the oscillating barrier, one can 
manipulate the stability of different states in a quantum potential. If in the static well there exists a bound state, then 
it is possible to stabilize a metastable state by adiabatically changing the oscillating frequency and amplitude of the 
barrier so that the unstable state eventually crosses over to the stable bound state. In this connection we note here 
that a related problem was studied in |Tsj . where the authors consider trapping an electron in thermal equilibrium 
within an infinite wall and a barrier which grows in time. There the confinement of particle in the well is achieved 
through the increase in the barrier height, whereas in our case it is accomplished mainly through the tuning up of the 
frequency of the oscillating barrier. 

In p"?! our main concern was the effect of avoided crossing on the quantum metastability in a periodically oscillating 
potential. In this paper we shall present a detailed study of the effect of direct crossing on the stability of a metastable 
system, and of how the oscillating amplitude affects the nature of direct /avoided crossings. We consider a simple 
quantum system in which a particle is trapped in between an infinite wall and a rectangular barrier. Two versions 
of the system are discussed. In Model A we consider the potential with an oscillatory height, while in Model B we 
make the bottom of the well oscillate instead. As in our previous work, the Floquet approach is adopted here. In each 
model, an equation which relates the Floquet quasienergy of the driven system to the amplitude and frequency of the 
time-periodic potential is obtained. This equation is solved numerically and the stability of the system is investigated 
as the frequency and amplitude of the oscillating field vary. To our surprise, the two models, very different in nature, 
turn out to possess the same Floquet energy spectrum (and thus the same stability) . A discrete transform connecting 
the two models is found which explains the exactness of the Floquet spectrum. A physical understanding of the 
equivalence of the two models is also given based on the principle of gauge invariance. 

The organization of the paper is as follows. In Sect. II, we will first briefly describe a simple metastable model of 
a-decay. Sect. Ill then introduces a modified model to include a time-periodically oscillating barrier (Model A). We 
will also describe the Floquet theory which we adopt in our analysis. The effects of the time-periodically oscillating 
barrier on the quantum metastability of the system is studied numerically in Sect. IV. Model B, which has the bottom 
of the well oscillating, is described in Sect. V. Sect. VI and Sect. VII present, respectively, the discrete transform and 
the gauge transformation that connect the two models. Sect. VIII summarizes the paper. 



II. SIMPLE STATIC MODEL OF QUANTUM METASTABLE SYSTEM 



A simple one-dimensional model to describe the phenomena of a-decay is that of a particle trapped in between an 
infinite potential wall and a finite barrier . The barrier potential of such a metastable system can be expressed as 



oo , x < , 
V{x) = { V > , a<x<b , 

, < x < a , and x > b 



(1) 



The wave function of the system has the well-known form ^>(x,t) = ip(x)exp(—iEt/h), where ip{x) satisfies the 
time-independent Schrodinger equation, 



d 2 



2m dx 2 

The wave function can be expressed in each regions as 



V(x) 



ip(x) = Etp(x) 



(2) 



Asm(kx) + B cos(kx) , < x < a , k=y2mE/h , 

= 4 Ce qx + De~ qx , a < x < b , q = ^2m(V - E)/h 2 

Fe lkx + Ge~ lkx , b < x . 



(3) 



Besides the ordinary boundary conditions which require that the wave function and its first derivative be continuous 
at the boundaries x — a and b, we could also impose another two boundary conditions on the wave function by physical 
consideration. One is B = as the wave function must vanish at the origin {x = 0); the other is G — since there 
is no extra nonvanishing potential in the region x > b to reflect wave from the right-hand side of the barrier. This 
later condition was first introduced by Gamow j2jj, Curney and Condon |2l|, and is usually called the "Gamow's 
outgoing wave boundary condition" in the literature. With these boundary conditions, we can eliminate the remaining 
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coefficients (A, C,D,F) of the wave function and obtain an equation that determines the eigenenergy of the system, 

% tan ka + 1 = f £ tan fca - l) e" 2 ^"") . (4) 

k q — ik \k ) 

By solving the above equation numerically for a given barrier height Vo, one finds that the solution of Eq.Q gives 
complex eigenenergy, E = Eq —iT/2. The imaginary part is related to the decay rate of the state. Despite some 
conceptual difficulties concerning the nature of complex energy, this simple approach by imposin g G amow's outgoing 
wave boundary condition successfully describes the exponential decay of metastable systems [l9L \22\ . 



III. MODEL A: THE STATIC MODEL WITH A TIME-PERIODICALLY OSCILLATING BARRIER 

We now modify the static model to include a time-periodically oscillating barrier, 

oo , x < , 

V(x, i) = { Vo + Vx cos(ut) , a < x < b , (5) 
, < x < a , and x > b , 

where Vo > and < V\ < Vq. Unlike the static case, in this time-dependent model we have to solve the time- 
dependent Schrodinger equation, 



ih 



dt 



2m dx 2 



*(*,*) ■ (6) 



A similar time-periodically oscillating potential (without the infinite wall) has been used to investigate electrons 
scattering in a time-dependent potential by Li and Reichl Q. They employed the Floquet formalism to set up the 
S-matrix of the driven system. The advantage of the Floquet formalism is that one could do away with any restriction 
of the amplitude and frequency of the driving force. Here we will mainly follow the procedure developed in to 
construct the wave function of our model of metastable system. 

The Floquet theorem states that the wave function of a time-periodic system has the form ^(x, t) = e _let ^$ e (x, t), 
where $ e (cc, i) is a periodic function $ e (a;, t) = & £ (x, t+T) with the period T = 2ir/u>, and e is the Floquet quasienergy, 
or simply the Floquet energy. It should be noted that the Floquet energy is determined only modulo Huj. For if {s, $ e } 
is a solution of the Schrodinger equation, then {V = e + nhui, <Ev = ®s ex-pjinuit)} is also a solution for any integer n. 
But they are physically equivalent as the total wave function ^> e is the same [23. All physically inequivalent states 
can be characterized by their reduced Floquet energies in a zone with a width Tiuj. We therefore consider solutions of 
e only in the first Floquet zone, i.e. e € [0, Tvjj). 

The Schrodinger equation inside the region of the oscillating barrier a < x < b is 



3 ft? 2 

ift— $ s (x, t) - — — $ e (x, t) + [Vo + Vi cos(cut)} * B (x, t) = e<S> e {x, t) . (7) 



Letting $ e (x,t) = g(x)f(t), Eq.Q can be separated into 

ft 2 d 2 g(x) 



2m dx 2 

and 



Vog(x) = Eg(x) , (8) 



iftffl _ Vl cos(wt)/(t) = (E- e)f(t) , (9) 
dt 

where E is a constant to be determined from the periodicity of the function f(t), i.e. f(t + T) = f(t). Eq.© can be 
integrated to give 

f{t) = e - i ( E - e ) t / h e^ iVl sin i ut )/ hw } (10) 

where we have set the initial condition as /(0) = 1. The second factor in Ea. ljlOp can be rewritten using the identity 


00 

e -iV 1S m{ut)/nu, = J n (a)e- mu]t ; a = -i , (11) 

^ — 4 Tiuj 
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where J n (a) is the Bessel function of the first kind. Ea. l|10|) then reduces to 

oc 

f(t) = e -^ E - e)t ' h J n (") e~ mut ■ (12) 

n=— oo 

The periodic condition f(t + T) — /(t), which reads 

oo 

f(t + T) = e -i(E-e)(t+T)/n J- J n (a) e - m ^ t+T ~> , (13) 

71 — — OQ 

OO 

= e -*(B-6)t/ft g J n ( a ) e -« e -i(E- e+n MT/R > ( 14 ) 
n— — oo 

= /(*) , (15) 

requires that (-E — e + nhu>)T/h = 2mK (m =integer). So the constant E can only assume the following values 

Ei = e + ZTiu; , I = m — n = integer . (16) 

The allowed values of E indicate that the time-periodic potential will generate side-band structures with Ei = e + ltkJ 
as the "side-band energy" of the l-th side-band. 

Knowing the allowed values of the constant Ei, we can now solve Eq.JHJ to get 



g(x)^J2 aieqiX + bie ^ lX > ?i = V 2m (^-£-N/?» 2 . (17) 
i 

Combining Ea. (|12[l and Ea. (|17|l . we obtain the wave function inside the region of the oscillating barrier a < x <b, 

n I 

= EE ( aie9ix + h ^ qix ) J «-< («) e ~ lEnt/h ■ ( 18 ) 



n / 



Having gotten the wave function inside the region a < x < b, the continuity of the wave function and its first derivative 
at the boundaries x = a and b at any time requires that outside the region a < x < b there are nonvanishing wave 
functions on each side-band channel. Therefore the wave function in the region < x < a is 



t) = J2 A n sm{k n x)e~ tE ^ h , k n = yf 2mE n /h 2 , (19) 

n 

which vanishes at the origin, ^/(0,i) = 0. On the other hand, the Gamow's outgoing wave boundary condition 
suggests that the wave function in the region x > b to have the form 

* m (x, f) = E t n e iknX e- iE ^ H . (20) 

n 

We can obtain relations among the coefficients A n ,a n ,b n and t n by matching the boundary conditions at x = a 
and b for each side-band. The results are 

A n sm(k n a) = ^ (a ie * a + b ie ' q ' a ) J„_,(o) , (21) 
i 

k n A n cos(k n a) = J2 91 (aie qia - bie- qia ) J n -i{a) , (22) 
i 

t n e tk " b = ^(a ie * b + 6 ; e-* b ) J„_,(a) , (23) 
i 

ik n t n e iknb = E * ( a ' e?ib " b ^ 9lb ) J «-'(«) • ( 24 ) 
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The Floquet energy is determined from these relations by demanding nontrivial solutions of the coefficients. In 
practice, however, we must truncate the above equations to a finite number of side-bands, eg., n = 0, ±1, . . -j±N . 
The number N is determined by the frequency and the strength of the oscillation according to N > Vi/Two 

To determine the Floquet energy e, we first separate the boundary conditions for the central band n = from the 
other side-bands I ^ in Ea. H21() - (|24|l (we set Aq = 1 for simplicity): 

sin(fc a) = (a e qoa + b Q e^° a ) J (a) + £ (cue*" 1 + he'*" 1 ) J_,(a) , (25) 
k cos(fc a) = q (a e* a - b e~«° a ) J (a) + £ qi (a,e^ - he^) J_;(a) , (26) 
i e*> b = (a e qob + boe^) J (a) + ^ (a;e*< 6 + 6 ie -* b ) J_ z (a) , (27) 

ifc * e ifc ° b = <z (a e*> 6 - & e-* b ) J (a) + £ « ( a ' e * b ~ 6 < e ^ ib ) J -'( a ) ■ ( 28 ) 

z#o 

The coefficients a z and 6; (Z ^ 0) in Ea. (|25l28|l can always be related to the coefficients ao and bo as 

a-i = fia(ko,w, Vi)a + fib(k ,u},Vi)b , (29) 
bi = gi a (k ,u),Vi)a + gi b (k ,u),Vi)b i (30) 

where /'s and g's are functions determined as follows. Eliminating the A n 's and t n 's in Ea. (|21|l - (|24|l . we can obtain 
A- n e«" a J a n + A+ n e-"" a J b n + £ ^,e*« J n _,o, + £ A+ je -" J n _ z b, 

i^n.O Z^n,0 

= -A- fi e«° a J n a Q - Al e- q ° a J n b , (31) 



and 



B+, n e^ b J Q a n + B- n e-i" b J b n + £ B^ h J n _ x a x + £ B~ l e~ Qlb J n _ib t 

-B+ e^ b J n a - B- Q e-*> b J n b , (32) 



where 



A^ l = cos k n a ± y- sinfc n a , and B^ = 1 ± i--^- . (33) 

/Crj, nJ^i 



As mentioned before, in practice a truncated version of Eqs. I|31|) and i|32|l has to be used. In a m side-bands 
approximation, the index n ranges from —m to m except 0. Then Eqs. (|31|l and (|32|l together represent a set of 2m 
equations with 2m unknowns a n and b n . The solutions of these unknowns can then be solved by the Cramer's rule 
Hg in the forms of Eq.® and Eci. (|5Ujl . 

Replacing the a;'s and 6;'s by means of Eas. (|29|) and (|3U[) . one can rewrite Eas. l|25|) to l|28|) as 

sin(fc a) = F 1 (ko;u,Vi)eP a a + F a (ko]u,Vi)e-*> a bo , (34) 

/c cos(/c a) - F3(*to;w,Vi)gbe* B ao-Jk(*o;w,Vi)«)e-* B 6o, (35) 

t e lkob = Fr^u^V^ao + Feiko^^e-^bo , (36) 

iMoe ifco6 = F 7 (fc ; w, V 1 )q Q e' lab a Q - F s {k ;uj, V^q^X , (37) 

where the coefficient Fi(ko] oj, Vi) (i — 1, ... ,8) are 

= Jo(a)+e-® a ^(/^ Q +Szae-* Q ) J_ ; ( a ) , (38) 
z#o 

F 2 (fc ;w,Vi) = J («)+ e 9oa ^(/ ;6 e 9ia +.gz^ a ) J_ ; (a) , (39) 

z#o 

F 3 (fco;w,^) = J ( a ) + e -^V « (/ iQ e^«- 5iQ e-« a ) J_,( a ) , (40) 
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F 4 (fc ;^^) = J ( a )- e ^V ^ {} lb e^ a - g lb e-^ a ) J_ ; (a) , (41) 
F B (*o;w,Vi) = JbW + e-^^^iae^+ftae-* 6 ) J_,(a) , (42) 
F 6 (fc ; W ,^) = J (a)+e 9ob ^(/^ b + g;b e-« fo ) J_,(a) , (43) 



2#o 



F 7 (fc ;w,^i) = Jo(a)+e^ nb V^(//ae 9lb - ffia e- 9ib ) J_ z (a) , (44) 
F 8 (fc ; W ,^i) = J (a)-e« o6 V^-(/ i6 e« ib - 9i6 e-«' fe ) J_;(a) . (45) 

Then by eliminating the coefhcients ao, bo and to in Eqs.(PU to l|37|l . we arrive at an equation which is a function of 

e (= h 2 k 2 /2m), 

F i —t&nk a + F 2 = — — — F 3 — tanfc a-Fi e S(H J . (46) 

ko F 7 q - iF 5 k \ k J 

Its solution determines the Floquet energy of the metastable system. Comparing this equation with the corresponding 
one for the static case, i.e. Eq.@, we see the effect of the time-periodic potential is contained entirely in the coefficients 
Fi. 

In the next section we will solve Ea. H46|) numerically. Once the Floquet energy e of the metastable system is 
determined, the probability of the particle still being trapped by the potential barrier at time t > 0, i.e. the non- 
decay probability P(t), can be obtained as 

= f Q b Mx,t)\ 2 d X 
f b \V(x,0)\ 2 dx 

= e 2Im(e)t/n fl \®s{x,t)\ 2 dx 

/>,(*, 0)| 2 d* 

with P(0) = 1. The imaginary part of the Floquet energy, which enters P(t) via the factor exp(2Im(e)t /%) , gives a 
measure of the stability of the system. Unlike the static case, however, here P(t) is not a monotonic function of time, 
owing to the time-dependent function h(t). But since h(t) is only a periodic function oscillating between two values 
which are of order one, the essential behavior of P(t) at large times is still mainly governed by the exponential factor. 
Hence, as a useful measure of the nondecay rate of the particle in the well, we can use a coarse-grained nondecay 
probability P(t) defined as |l7| 

P(t) = e 2Im{€)t/h {h(t)) , (48) 

where (hit)) is the time average of h(t) over one period of oscillation. 

Owing to the complexity of the coefficients F^ko; ui, Vi), the solution of Ea. (|46|l has to be obtained numerically. We 
shall present the numerical analysis in the next section. Here let us consider the limiting cases of our model. From 
the expressions of the coefficients Fi(ko] oj ,V\) fEa. (|38fl to Ea. (|45[) ). we can easily check the coefficients Fi(ko;u>, Vi) 
all approach one when the parameter a = Vi/Tiu — > 0, 

lim F,(fc ;w,Ui) — > l,i = l,..., 8. (49) 

Vi/ftu;->0 

Hence in the limit V\ — ► or ui — > oo, Ea. (|46l) is identical to Eq.@ in the static case, and the Floquet energy in 
these limits are just the complex eigenenergy of the static case. This is understandable, since in the limit V\ — ► the 
potential becomes static. And at very high frequencies, the dynamical time scale of the particle in the well is much 
larger than the time scale of the oscillating barrier hence the height of the time-periodically oscillating barrier 
seen by the particle should be its time-averaged value, namely, < Vq + V\ coswi >=Vq, which is independent of the 
time-oscillating part of the barrier [j| . 
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IV. NUMERICAL ANALYSIS OF MODEL A 

In this section, we will investigate the quantum metastability of Model A by finding the solutions of Ea. l|46|) 
numerically. We are interested in the metastability of the system under different amplitudes and frequencies of the 
oscillating barrier, since these are the experimentally controllable parameters. For definiteness, we take the barrier 
to have height Vq = 10.0 in the atomic units (a.u.) (e = m e = % = 1), and boundaries at a = 1 and b = 2 (a.u.). 
With these parameters, there is only one solution with Re(E) < Vq in the static case. This is a metastable state 
having complex energy Eq/Vo = 0.322052 — 0.000110412i. Other solutions have Re(E) > Vq and are more unstable. 
We shall be interested in this paper only in the coupling between the metastable state Eq and the next higher one 
with -Ei/Vo = 1.11205 — 0.025062i. Coupling between other metastable states can be considered similarly. For 
the oscillating potential, we shall solve Ea. (|4*H)l numerically for cases with different oscillating amplitude Vi/Vq and 
frequency lu/Vq. 

A. Comparison of 2 and 3 side-bands approximations 

In this subsection we first study how good a 2 side-bands approximation is by comparing it with the 3 side-bands 
approximation. Floquet energy as a function of the amplitude of the oscillating barrier V\ evaluated based on 2 
side-bands (solid line) and 3 side-bands (dotted line) approximations are shown in Fig.(l) and Fig. (2) for lo/Vq = 0.01 
and 0.02, respectively. 

As mentioned before, if {e, $ e } is a solution of the Schrodinger equation, then {e 1 — e + rihuj, $ e ' = $ £ exp(inujt)} 
is also a solution for any integer n. And they are physically equivalent as the total wave function is the same. For 
metastable system, this means that all Floquet states with real parts differing by nhui (n — 0, ±1, . . .) will have the 
same imaginary part. We have checked this for n = 0, ±1 in the 2 and 3 side-bands approximations. 

Part(a) [part(b)] of Fig. 1 and 2 show the relation of the real (imaginary) part of the Floquet energy and the 
oscillating amplitude parameter Vi/Vo at fixed ui. Fig. (1-a) and (2-a) show that if Re(e/V~o) = Re(E n /Vo) = 0.322052 
is the Floquet energy of the system, so are the Re(e/Vo) ±ui/Vq for all values of V\ less than certain limit. Fig. (1-b) 
and (2-b) indicate that the imaginary part of these solutions merge into one curve for V\/Vq < 0.01 in Fig. (1-b) and 
Vi/Vq < 0.02 in Fig. (2-b). The divergence of the three curves at higher V\ is caused by the truncation of the number 
of side-bands in numerical analysis. Fig. (1-b) shows that the 2 side-bands approximation is good for V\/Vq < 0.01 at 
ui/Vq = 0.01 (i.e. a = V\/w < 1), and the 3 side-bands approximation is accurate for Vi/V < 0.02 at ui/Vq = 0.02 (i.e. 
a = V\/uj < 2). These results are consistent with the usual requirement for the N side-bands approximation mentioned 
previously, namely, a = Vi/lu < N [71 la. uj\. Similar argument for the accuracy in the approximation can also be 
checked in Fig. (2-b). Stability of the system is determined by the imaginary part of the Floquet energy. From part (b) 
of Fig. 1 and 2, we see that the imaginary part of Floquet energy increases as V\ increases. This means that, at fixed 
frequency, a larger oscillating amplitude will cause the system to be less stable. On the other hand, once V\ is reduced 
to zero the quantum metastability of the system will approach to the static case {Log(— Im(e/Vo)) = —3.95698), 
because, as we explained at the end of the last section, in the limit V\ — * the driven system reduces to the static 
case. 



B. Floquet energy as function of the frequency and amplitude of the oscillating field 

In this subsection, we will investigate the effects of the oscillating frequency and amplitude on the quantum metasta- 
bility of the system in the 2 side-bands approximation. As we have learned in the previous subsection, the 2 side-bands 
approximation is good enough for uj > V\ (a -1 > 1). 

Fig. (3-a) and (3-b) present, respectively, the relation of the real and imaginary parts of the Floquet energy and 
the driving frequency with the parameters Vq = 10.0, a = 1, b = 2, and V\ = 1.0 in the atomic units. The solutions of 
Ea. (|46|l have the form e = Eq + nuj (n = 0, ±1, ±2, . . .), with Re(eo) (the horizontal branch) lie close to the energies 
Re(Eo) and Re(Ei) in the static potential. That is, these branches of Re{e) emanate from either Re(Eo) or Re(Ei) at 
uj = 0. Branches emerging from the same point have the same imaginary part. Numerical results show that, with the 
barrier oscillating, these states become less stable. For simplicity, in Fig. (3a) we show only two branches emerging 
from the two states with Re(Eo) (solid curves for n = 0, 1) and Re(Ei) (dotted curves for n = 0, —1) in the static 
case. As mentioned before, solutions can always be reduced to the first Floquet zone, Re(e) (modulo u>), which are 
points under the line Re(e) = uj. 

Figure (3-b) shows that the imaginary parts the Floquet energies of the two states fluctuate slightly as the frequency 
changes. At frequency close to the resonance frequency uj pa lo r = Re(Ei) — Re(Eo) w 0.79Vb, a direct crossing of 
the real parts of the Floquet energies of the two states occurs. This is indicated by the point a in the first Floquet 
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zone, or equivalently, point b in the second Floquet zone. That it is a direct crossing is shown more clearly with a 
refined scale in Fig. (3-c) and (3-d). We see that for the first metastable state (solid curve), Re{e) exhibits around 
lor a "Fano" resonance pattern, where a sharp dip is followed by a peak. For the second metastable state (dotted 
curve), an "anti-Fano" pattern was seen instead. In the neighborhood of ojr, the first metastable state becomes less 
unstable (point c), while the second state appears to be slightly more stable (point d). Beyond lur, the stability of 
the two states revert to their respective orders of magnitude. Similar behavior is also observed at the point of direct 
crossing in the system discussed in |17|. 

Figures 4 to 7 present similar graphs as in Fig. 3, but with Vx/Vq — 1.6,1.7,2.0 and 4.0, respectively (other 
parameters being the same). In Fig. 4 we see that, with Vi/Vq increases up to 1.6, the system still behaves essentially 
in the same way as when Vq/Vi = 1.0. Only the first metastable state becomes less stable, while the second metastable 
state tends to be more stable at the point of direct crossing. 

The behavior of the system changes qualitatively, when Vo/Vi increases beyond a critical point. At Vi/Vq = 1.7, 
the crossing between the two states changes from a direct crossing into an avoided crossing, as depicted in Fig. 5. 
Increasing Vi/Vq further only enhances the difference of the Floquet energies at the point of avoided crossing (Fig. 6 
and 7). At and beyond the avoided crossing, the two states exchange stability: the second and originally less stable 
state becomes more stable than the first state, which now assumes a stability which is of similar order of magnitude 
as that of the second state before the crossing. Such phenomenon has been previously discussed in |l7j . 

The above results indicate the possibility of controlling stability of states in a quantum well. One may tune up 
the frequency adiabatically up to the point where the real parts of the Floquet energy of the two states meet. If 
the crossing is a direct crossing, then one can increase the magnitude of the oscillation V\ of the barrier before the 
crossing until the direct crossing turns into an avoided one. When the frequency is increased further to pass beyond 
the avoided crossing, the two states exchange stability. Then by reducing V\ to zero, the two states in the static well 
are interchanged. If the lowest energy in the static well is a stable bound state (as in the case discussed in ^3)> the 
above procedure can turn an upper metastable state into a stable one. 



V. MODEL B: THE STATIC MODEL WITH A TIME-PERIODICALLY OSCILLATING BOTTOM 

Let us now consider a variant of the above model, namely, the static model with a time-periodically oscillating 
bottom instead of an oscillating barrier. The potential can be expressed as 

{oo , x < , 

£cos(arf), 0J*<a, (50) 
Vo , a < x < b , v ' 

, x > b . 

The wave function inside the barrier well, i.e., < x < a, can be obtained by the similar procedure as in Sect. Ill, 
except that we have to impose the boundary condition $f'r(x = 0, t) = in this case. The result is 

%{x,t) = ^^A; i sin(fc^)J i _ n (a)e- J ( £+m ^ t / ?i , 

l n 

= J2Y. A 'n s ™(knx)Ji- n (a)e- lE ^ h , (51) 

l n 

where k n = y^2m(e + nhu>)/h. Continuity of the wave function at the boundaries x — a and b requires that the wave 
function has the form 

*//(M) = ( a i e9 ' X + h 'i e ~ qiX ) e~ lElt/h (52) 
i 

in the region a < x < b, where qi — y^2m(Vo — e — lhcu)/h, and 

*izr(*> *) = tle^e-^'^ , (53) 
i 

in the region x > b, where the Gamow's outgoing wave boundary condition has been imposed. 

The relations among the coefficients A' n ,a' nl b' n and t' n can be obtained by matching the boundary conditions at 
x = a and b. This leads to 

J2<Mk n a)Ji- n (a) = a[e^ a + b[e-^ a , (54) 
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Y,kiA' n cos(k n a)J^ n (a) = m (a'^' a - b'^) , 

n 

t[e iklb = a 1 ^ + b' l e- qib , 
ifc^e**' 6 = qi (ale^ - b[e-^ b ) , 

Eliminating t\ from Ea. (|56fl and Ea. (|57fl gives the relation between a[ and b\, 

= _l + lqi / kl 
1 1 - iqi/k 1 

Using this relation to replace b[ in Eq. l|54[l and Ea . (|55|) . we obtain 

~ V l-iqi/h J 

and 



A' n k n cos(k n a)Ji- n (a) = qie q 



lib ( - qi {b- a ) , l + iqi/h qi{b - a) \ , 



I ■ 



1 - 

respectively. These two equations can be combined by eliminating the coefficients a',, giving 



(55) 

(56) 
(57) 



(58) 



(59) 



(60) 



E 



B + 

A - e -«(6-a) M A + ?,(6-a) 

n,t t->— n,i 



k n Ji- n (a)A' n = 



(61) 



where the notation ^4^ z and are as defined in Model A, Eq. 

We now show how to determine the Floquet energy of the metastable system from Eq. I|61|l . For a general amplitude 
Vi 7^ 0, Ea. (|61fl in the m side-bands approximation is 



E 



B + 

A - e -qi(b-a) U A + qilb-a) 



B 



1,1 



k n Ji- n (a)A' n = 



(62) 



Moving the term with the coefficient A' of the above equation to the right-hand side, we get 



E m 

n#0 



B + 

A- p-1l(b- a ) hi A + p <ll{b-a) 



k n Ji-n(u)A' n 



B 



A - -mib-a) _ _J£ A + qi (b-a) 



k Ji{a)A' , 



(63) 



The left-hand side of the equation contains 2m unknown coefficients A' n ^ . As the total number of side-band is also 
2m (I = ±1, ■ • ■ , ±m), it is sufficient to use the Cramer's rule to obtain the relation 

^ = C n (e,7i,«K ( 64 ) 
with Co = 1. We now insert this representation of A' n _^ into Ea. (|62[l . Ea. l|62|) for the central band (I — 0) becomes 



E 



A- -qo(b-a) °-° A+ p qo(b-a) 



B, 



0.0 



k n J- n (a)C n } A' = . 



(65) 



Nontrivial solution of A' , i.e. A' ^ 0, requires that 



E 



n——m 



B + 

A- p-qo{b-a) °'° 4+ q (b-a) 
-E>n n 



k n J- n {a)C n = 



(66) 



The Floquet energy of the metastable system e can then be determined by finding the solutions of this equation. 
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Let us mention briefly the limiting case in which Vi — ► (i.e., a — > 0). Since Jo(a = 0) = 1 and Ji^o{a = 0) = 0, 
Eq.JSSJ) reduces to 

B + 

A o e~ qo{b ~ a) - -^-A+ e q ^ b - a ^ = , (67) 
B o,o ' 

which, as we expect, is just the result of the static case 

( Si tan k a l) e^-) = |° tan/fc a + 1 . (68) 
q - ik \k Q J k 

We had solved Eq. (|dd) numerically for the Floquet energy of the system. To our surprise, the functional dependence 
of the Floqeut spectrum on the amplitude and the frequency of the oscillation turns out to be exactly the same as 
that of Model A. Since the sets of equations of boundary conditions determining the Floquet energy are so different in 
form for the two models , that they should give exactly the same solutions can be of no coincidence. This prompted 
us to search for the underlying connection between the two models. After some efforts we eventually realize that these 
models are related by a discrete transform in the sense that Eqs. ((21|) - (I24() can be transformed into Eqs. H54|) - (|57|l . and 
vice versa, by the discrete transform to be described below. Since Eq. (1221) [ 1)55(1 ] and Eq. ((24(1 [ 157(1 ] are obtainable 
from Eq. 1(211) [ l(54(l ] and Eq. 1(231) [ l(56(l ] by differentiation with respect to the parameters a and b, respectively, it suffices 
to show the transformations between Eqs. ((21() and 1(54(1 . and between Eqs. 1(23(1 and ((56(1 . This will be discussed in 
the next section. More recently, we have also realized that the underlying physics responsible for the equivalence of 
the two models is the gauge equivalence of the models, which we shall present in Sect VII. 



VI. DISCRETE TRANSFORM 



The aforementioned discrete transform is defined as follows. Consider an infinite discrete sequence {g n } = 
{. . . , <7_i,<7o> git ■ • •} of complex numbers g n (n = — oo, . . . , oo) with finite norms. Let us form another sequence 
{g' n } by the following transform based on the Bessel functions Jfe(a) with integral order k and argument a: 



9i («) := Hi [{gn},a\ 

OO 

= Yl (-V n 9nJi-n{a) , n = 0,±l,±2,. 



(69) 



This transform has the interesting property that its inverse transform is given by the same transform, i.e. the sequence 
{g n } is obtainable from {g' n } by the same transformation: 

g n = H n [{g'i\,a\ 

OO 

= J2 i-l) l 9iJn-i (a) ■ (70) 

I — — OO 

To prove this, let us substitute Ea. l|69|) into the right-hand side of Eq. J70J)- This gives 

oo 

]T (-i)%'J„_K«) 



l=—oo 



= E E (-l)"W*-m («) 

Z— — oo [_m— — OQ 

oo oo 

= J2 (- 1 ) m 9™ E (-!)' J '-™ ( a ) J »-< («) 

m— — oo / — — oo 

oo oo 

= E (-!) m+n 5™ E Jl-m{a)Ji-n{oc) ■ 



a) 



(71) 



/ — — OO 



In obtaining the last line in Eq. <f7T|> . we have made use of the property J n -l{oc) — (— Jl-n{pc). To proceed 

(72) 



further, we make use of a version of the addition theorems of the Bessel functions, namely [24 

oo 

^ Jk+v ipt\) J k (02) = Jv {pL\ - a 2 ) . 



k= — oo 
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It is easy to prove that this theorem leads to the following interesting identity for Bessel functions with integral orders: 

oo 

^ Jl-m («) Jl-n («) = Jm-n (0) 

I — — OO 

= $mn ■ (73) 

Using Eq. (|73l) . Eq. (|71|l becomes g n , which is exactly the left-hand side of Eq. (|70|l . 
Let us now apply the transform to Eq. I|54|) . We get 

(-l) n A' nS m(k n a) = £(-1)' (ale^ + ble-^) J n _,(a) . (74) 
i 

We recognize that this is the same equation, Eq. (|21[l. satisfied by j4 n , a„ and 6 n , with the same Floquet energy £ 
and frequency w. This implies that the two sets of coefficients are related by 

A n = (-l) n A' n , a n = (-l)X , b n = (-l) n b' n . (75) 

Similarly, by applying the transform to Eq. 121(1 . we get Eq. I|54|l with the same connections, namely, Eq. (|75|) ■ among 
the coefficients. With this connection of the coefficients it is seen that Eqs. (|21[1 and l|54|l are dual pairs under the H 
transform. 

Now consider the transform of Eq. (|56|l . With Eq. I|75|) , we get 

5>l)^e* b J_ ; (a) = (a^ + ^e-^) ./„_,(«), 

i / 

= ^(a^ + V^) J n _,(a) . (76) 

This is the same as Eq. t\2'M with the same Floquet energy and frequency, provided that the coefficients t n and t' n are 
related by the W-transform in the following forms: 

oo 

t n e^ b = (-l)'*je tt,6 J»-i (") - 

I— — OO 
OO 

= J- (-l)'i/e* b J„_/(a) . (77) 

/— — OO 

The above discussions show that, the two models are mapped into one another under the 7i-transform, with the 
relevant coefficients in the wave functions being related through Eqs. (|75|l and (|77|l . Hence the Floquet energy spectra 
as function of oscillation's amplitude and frequency of the two models are identical. The wave functions, however, 
are different, since they have different forms, namely, Eqs. l|18|) - (|2Uf) for Model A, and Eqs. H51 M 53 fl f° r Model B. 
Nevertheless, it turns out that the relations between their coefficients in the region < x < b, where the particle is 
confined, can be understood as the result of a gauge transformation of the two wave functions. This will be discussed 
in the next section. 



VII. GAUGE EQUIVALENCE OF THE TWO MODELS 

In this section we show the underlying connection between the two models based on gauge invariance principle. 

One needs only to consider the region in < x < b , since the metastable system in these two models is mainly 
confined within this region. The region x > b does not affect the essential physics underlying the two models. This is 
because for a metastable state, its wave function is extremely small (nearly identical to zero, as confirmed numerically) 
in the region x > b before the system has actually decayed. This region serves only as an ideal drain. 

The gauge transformation of a quantum system with an ^-independent gauge function \(t) is given by: 

ip(x,t) i>' = e^ x(t) V>(M) , 

V(x,t) -> V'(x,t) = V-%- X ® ■ (78) 
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When the gauge function x(t) is chosen to be 



X(t) = — sin(wt) , (79) 



the potential in region < x < b in Model A, Eq. (|5J), is transformed into 

nM) = (r/ lCOs( ^' °i x <V (so) 

Vo j a < x < b . v/ 

This potential is physically equivalent to the potential describing Model B in the same region: they only differ in 
the origin of the time variable t by an amount ir/uu. Hence, the potentials of the two models are related by a gauge 
transformation with gauge function Eq. I|79|l and a time shift t — > t + tt/u. 

We now examine the connection between the wave functions of the two models in the regions considered. Under 
the gauge transformation described by Eqs. 178|) and 179fl . the wave function of Model A transforms according to 

(x, *)-►#' = e lVl sinM)/^^ t j ^ ^ 

with 



iVi sin(wt)/ 



huJ = J2 {-^) m Jm{ot)e- imut (82) 



m— — oo 



which can be easily obtained from Eq. I|ll|) by changing V\ to — V% and using the fact J n (—a) = (—l) n J n (a). 
According to the previous discussion, we shall make the shift t — > t + tt/lu in and compare the resulted wave 
function W(x, t + tt/uj) with that of Model B. 

First we consider the region < x < a. The wave function of Model A in this region, Eq (jl9[) . is gauge-transformed 
into 

*J(a:,t) = e lVlsin ^ t ^ hL °^2 A ^ a[l, -(knx)e- lE " t/h 

n 

= E ^2(- 1 ) mj m (<*) An sin{k n x)e- l{E+{n+m)huj)t 

m n 

= E E(- 1 )^" J '-« ( a ) A ™ sm^K^ 4 ^ . (83) 

n I 

In the last line of Eq. 1(83(1 the index m is changed to m = I — n, and Ei — e + (m + n)fiuo. With the shift i — > t + 7t/cl>, 
the time-dependent factor in Eq. 1)83(1 picks up a factor exp(~iEit/h) = (— 1)' exp(— iire/Tiuj). Accordingly the wave 
function is shifted to 



*j(af,t) = V'(x,t + n/u) 

= e-***"*' J2 E^ 1 )"^" sin(*n«)Ji-» («) e" lB,t/?i ■ (84) 



i n 



It is now clear that the wave function 'J" in Eq. (|84|l is equivalent, up to an irrelevant phase factor, to ^>'j of Model 
B, Eq. I|51(l . with the identification of A' n = (—l) n A n as discovered in the last section. 

Now we turn to the wave functions in region a < x < b. The relevant wave function of Model A, Eq. ((18(1 . is 
gauge-transformed to 



-i(E n +mKw)t/H 

E E E ( a ' e ** + 6 ' e "^) J ™-' ( a ) J »- fe ( a ) e" ,;Bfc * /?l , (85) 



k n I 



where we have used J- m ( a ) = ( — ^-) m Jm (oc), and relabeled the index in the last expression. By making the same 
shift in the time variable, the wave function changes to 

= e-/'- ^ ^(-i)^ ( a , e ^ + 5 ie -^) e -* f / s E 4-i H J n -k(a) . (86) 

k I n 
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The sum over n in the last factor is simply the Kronecker delta symbol Sik according to Eq. (|73(l . Hence 

V'jjix, t) = e -^/^J2(-l) 1 (aie qix + he- qiX ) e ~ lE ' t/h . (87) 
l 

Again, up to the irrelevant phase factor, the wave function ^"j is equivalent to of Model B, with aj = (— 1) a; 
and b[ = (—l) l bi as before. 

The gauge equivalence of the wave functions of the two Models in the physically relevant regions provides a physical 
understanding of the identity of their Floquet spectra. 

Finally, we note here that from such gauge invariance argument it also follows that the static problem is equivalent 
to the problem with the bottom of the well and the top of the barrier oscillating in phase with the same frequency 
and amplitude. 

VIII. SUMMARY 

In this paper we investigate the stability of a particle trapped in between an infinite wall and square barrier, with 
either a time-periodically oscillating barrier or bottom of the well. In these models we do not restrict the amplitude 
of the time-periodically oscillating potential to be small, so a nonperturbative approach based on the Floquet theory 
is adopted to tackle the problem. For each case we derive an equation to determine the Floquet energy of the system. 
The imaginary part of the Floquet energy can be used as a measure of the stability of the system. It is found that in 
general the oscillating field causes a quantum metastable state to decay faster. But when the Floquet energies (real 
part) of two states make a direct crossing, the more stable of the two states will undergo a resonance-enhanced decay, 
whereas the less stable state is slightly stabilized near the resonance frequency. Beyond the resonance frequency, the 
states resume stability of the original order of magnitude. However, if the crossing is an avoided one, the two states 
will switch stability: the less (more) stable state becomes more (less) stable. This gives the possibility to stabilize a 
metastable state, as discussed in |17| . Here we further found that when the amplitude of the oscillating potential is 
increased, a direct crossing could change into an avoided one. Thus, our results show that one can manipulate the 
stability of different states in a quantum potential by a combination of adiabatic changes of the frequency and the 
amplitude of the oscillating barrier. The two models we considered here were found to have identical Floquet energy 
spectrum for the same frequency and amplitude of the oscillating field. This was explained by a discrete transform 
which connects the equations of boundary conditions of the two models, and by a gauge transformation which maps 
the wave functions of the two models. These arguments also guarantee that the model discussed in [TtJ will have the 
same Floquet spectrum, whether the oscillating part is the barrier or the bottom of the well. 
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Figure Captions 

Figure 1. Relation of the real (a) and imaginary (b) parts of the Floquet energy vs. the parameter V\/Vq at fixed 
frequency parameter lj/Vq = 0.01 for Vo = 10.0, a = 1.0,6 = 2.0 (a.u.). Solid and dotted lines represent the 
values evaluated in the 2 side-bands and 3 side-bands approximation, respectively. Values of Vi/Vb in (b) beyond 
which the lines diverge indicate the points of breakdown of the respective approximations. 

Figure 2. Same plots as Fig. 1 for Vo = 10.0, a = 1.0, b = 2.0 (a.u.) and the fixed frequency parameter uj/Vq — 0.02. 

Figure 3. Relation of the real and imaginary parts of the Floquet energy vs. the the driving frequency oj/Vq for 
Vo = 10.0, V\ = 1.0, a = 1.0 and b = 2.0 (a.u.). Fig. (a) and (b) show the real and imaginary parts, respectively, 
of the Floquet energy. Fig. (c) and (d) depict the real part of the Floquet energy of the lower and the upper 
state, respectively, at the direct crossing with a refined scale. 

Figure 4. Same plots as Fig. 3, but with V\ = 1.6 (a.u.). 

Figure 5. Same plots as Fig. 3 with V\ = 1.7 (a.u.). Note that the crossing has become an avoided one. 

Figure 6. Same plots as Fig. 5 with V\ = 2.0 (a.u.). 

Figure 7. Same plots as Fig. 3, but with V\ = 4.0 (a.u.). The avoided crossing is enhanced. 
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